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The dependence of the size of the thermal source on the centrality in 
ultrarelativistic heavy-ion collisions is studied. The interaction region con- 
sists of a well defined thermalized core, and of an outer mantle where the 
production scales with the number of participants. The thermal source 
builds up in the region with the largest density of participants in the trans- 
verse plane. Particle production in the thermalized core is enhanced in 
comparison to the wounded nucleon model. The change of the degree of 
strangeness saturation with centrality is also discussed. We perform an es- 
timate of high p± jet absorption finding that an increase of the absorption 
in the thermal core is compatible with the data. 

PACS numbers: 25.75.-q; 24.85. +p 

1. Introduction 

Relativistic heavy-ion collision experiments at RHIC energies y/s = 
130GeV and 200GeV reveal a large degree of thermalization in the produced 
system. Particle spectra have thermal shapes, with a substantial amount 
of collective flow. The presence of a strong transverse and elliptic collec- 
tive flow demonstrates the creation of a strongly interacting phase in the 
history of the system. This stage can be described to a large extent by an 
ideal hydrodynamical evolution. Particle production in heavy-ion collisions 
is characterized by the thermal parameters of the source : its tempera- 
ture, transverse flow and possibly the degree of strangeness saturation. The 
temperature, the chemical potentials and the strangeness saturation at the 
chemical freeze-out can be extracted by fits of a statistical hadronization 
model to the ratios of abundances of particles produced in the central ra- 
pidity region. The description of particle spectra in transverse momentum 
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requires the knowledge of the collective flow profile in the source and of the 
temperature at the kinetic freeze-out. 

Usually it is assumed that all the produced particles are emitted from a 
single thermal source. One finds that the parameters of the thermal freeze- 
out depend on the centrality. For the chemical freeze-out one observes a 
significant variation of the strangeness saturation factor in peripheral colli- 
sions p] 12] . For the kinetic freeze-out a stronger collective flow and a lower 
freeze-out temperature comes out from the fits to the spectra of particles 
in central collisions A different approach allows for spatial variations 
of the parameters of the freeze-out surface surface, i.e. temperature and 
chemical potential gradients [I]. 

The extraction of the thermal parameters of the source and the descrip- 
tion of its dynamics in hydro dynamical models are useful in the description 
of the evolution of the dense matter produced in the collision. Such studies 
include the modeling of the elliptic flow in non-central collisions. Quantita- 
tive results require the knowledge of the initial conditions in the simulation. 
Therefore the understanding of the variation of the size and of the nature 
of the thermal source with centrality is crucial. 

2. Size of the thermal source 

Experimental data at the highest RHIC energies demonstrate the forma- 
tion of a thermalized system in Au + Au collisions : the spectra are thermal 
with a collective transverse flow component, the elliptic flow is compatible 
with an ideal fluid dynamics, particle ratios can be described by universal 
chemical freeze-out parameters. From the theoretical side we expect that at 
high energy densities formed in the collision the system undergoes a phase 
transition to the quark-gluon plasma. This phase transition could facilitate 
a complete chemical equilibration. If the partonic plasma is strongly inter- 
acting near the critical temperature, the system would behave as an ideal 
fluid. These feature should not vary with centrality if the thermalization 
and the evolution of the thermal source are dominated by the phase transi- 
tion. The kinetic freeze-out is given by the local density and flow gradient 
and should not vary strongly with the impact parameter of the collision 
either. 

Therefore a reasonable assumption would be to have a thermal source 
with the same thermal parameters at all centralities; the variation of the 
particle production coming only from the change of the size of the source. 
It has been noticed that particle production per participant in heavy-ion 
collisions at RHIC energies is significantly larger than in p + p or d + Au 
collisions [SI [Zl [HI - m the analysis of the experiments at RHIC energies 
Glauber Monte Carlo calculations based on the wounded nucleon picture 
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permit to extract the number of participants (or equivalently wounded nu- 
cleons) used to quantify the difference of the production per participant in 
Au + Au or d + Au collisions with respect to p + p interactions. Geometri- 
cal scaling predicts that the bulk of the particle production, dominated by 
soft particles, scales with the number of wounded nucleons [5]. It has been 
noticed that the density of charged particles at r/ = and the total number 
of charged particles produced in d + Au collisions scales with the number of 
participants : 

dN dAu dNp P 

~*r - Npart ^r ■ ( } 

Remarkably, this scaling is fulfilled in a more general way by the distri- 
butions of charged particles in a broad range of pseudorapidity 10 . On 
the other hand particle production in Au + Au collisions deviates from a 
simple superposition of independent production from all the wounded nu- 
cleons. The experimental data show a larger multiplicity than predicted by 
the scaling of the p+p production by N part /2. This and the evidence of the 
thermalization discussed above point to the importance of the second stage 
of the collision. The main characteristic of this stage is the existence of a 
thermally equilibrated source. Particle production from the thermal source 
occurs late and the production mechanism is far from the superposition of 
elementary collisions. In a system undergoing Bjorken longitudinal expan- 
sion the number of particles emitted at central rapidity is proportional to 
the transverse extension of the source. If the density of elementary col- 
lisions in the transverse plane is high one expect fast thermalization and 
eventually a phase transition to the quark-gluon plasma. In a geometrical 
Glauber model the decisive parameter is the density of wounded nucleons 
in the transverse plane. In peripheral collisions the density of colliding nu- 
cleons is small, and the production would be similar as in d + Au collisions. 
In semi-central and central collisions, at the center of the interaction re- 
gion the density of wounded nucleons is large and the transverse energy 
deposited is high enough to guaranty a fast thermalization. The boundary 
of the thermal source is given by the limiting density of wounded nucleons in 
the transverse plane sufficient for the fast thermalization processes to occur. 
The observation of universal chemical freeze-out parameters in Au + Au col- 
lisions and the indication of the formation of a strongly interacting plasma 
at near critical conditions, lead to the assumption that the boundary of the 
thermal source and of the quark-gluon plasma phase are the same. In the 
outer part of the interaction region the density of participants and therfore 
the density of the deposited energy is not high enough to guaranty the ap- 
proach of the critical point an the resulting fast thermalization. At RHIC 
energies the boundary in the interaction region between the thermalized 
source and the outer mantle is well defined by the criterion of the phase 
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transition. It turns out that the density of participants in the outer man- 
tle is low, they undergo at most 2.5 binary collisions per wounded nucleon 
pair. A similar number of binary collision is predicted for central d + Au 
collisions. The number of particles produced in d + Au interactions (even 
in the most central collisions) scales with the number of participants. The 
low number of binary collisions in the outer mantle justifies the assumption 
that the number of charged particles produced in this way is proportional 
to the number of wounded nucleons in the mantle N™ rt . 

The number of particles emitted from the source depends on the source 
size at the freeze-out. Since the freeze-out conditions are universal, the 
final number of particles produced should scale with the transverse energy 
deposited in the thermal source. The deposited energy is proportional to 
the number of wounded nucleons Np art in this high density core of the 
interaction region. The density of charged particles produced in the central 
pseudorapidity region is 

dNck _ N^ndN% + N^dN^ ^ 



drj 2 di] 2 drj 

•p 

dr\ 



dN pp 

where — is the density of particles produced produced in p + p colli- 



dN th 

sions and is the density per participant pair from the thermal source. 
The enhancement of the particle production in the thermal source can be 
effectively described by one parameter a 

The number of wounded nucleons in the mantle is calculated from the 
formula 

N ~ H = j fgfgzn^ _ • w 

where d ^2 rt is the density of participants in the transverse plane, calcu- 
lated from the Glauber Monte Carlo model. The integration is restricted 
to densities of participants lower than d cu t and involves only particles emit- 
ted in directions other than the dense core (the function @((f>) is one if the 
emission direction is not shadowed by the core and zero otherwise). The 
number of participants in the core is calculated analogously 

A£k = /Ag^ty*- -%*)(!-«»<#)) 

f 2 d(f> d 2 N part d 2 N par t , \ 

+ r s 2^^ ( ^p — d ™ th (5) 
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Fig. 1. Density of charged particles produced at central pseudorapidity in Au + Au 
collisions scaled by the density of particles in p + p collisions from the BRAHMS 
[TTIIH, PHENIX and PHOBOS |Sl|7j experiments compared to the predictions 
of the core-mantle model (solid line). The dashed line represents the number of 
participant pairs. 



the main contribution comes from the region with the density of participants 
larger than d cu t- For the calculation of the density of wounded nucleons we 
use Woods-Saxon nuclear density profiles 1/ (1 + exp((? — Ra)/cl)), with 
parameters a = 0.535fm and Ra = 1.12^4 1//3 fm. The cross sections are 
a = 42mb and 41mb for ^/s = 200GeV and 130GeV respectively. 

The results of the calculations for the scaled particle pseudorapidity 
density 

dNAuAu i dNpp 
drj drj 

are presented in Figure ^ Published experimental data from the BRAHMS, 
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PHENIX and PHOBOS collaborations are shown. The density of charged 
particles divided by the p + p density increases faster than the number 
of participant pairs shown by the dashed line. Moreover this increase is 
not exactly linear in N par t. There is a strong increase of dN/dr] around 
30 — 100 participants and then a gradual increase of the slope up to the 
highest centralities. We can reproduce this effect in the core-mantle model 
with the parameters d cu t = 2fm -2 and a = 0.65 at y/s = 200GeV and 
d C ut = 2.2fm -2 and a = 0.5 at yfs = 130GeV. At centralities corresponding 
to Np ar t = 40 the thermalized core appears in the interaction region, caus- 
ing a fast increase of the scaled particle density. For more central collisions 
the contribution of the thermal core in the total production increases grad- 
ually up to 95% at zero impact parameter. This leads to a slightly faster 
than linear increase of the particle density with the number of participants, 
similar as seen in the data. The increase of the production per participant 
pair in the core 1 + a and the decrease of the cutoff density of participants 
d C ut with energy indicate that the wounded nucleons get effectively "fatter" 
in the the transverse plane and deposit a larger transverse energy at higher 
collision energy. At SPS energies the mantle part of the interaction region 
is dominant. The number of binary collision in the mantle increases and 
the assumption of the N™ rt scaling of the production in the mantle breaks 
down. 



3. Strangeness 

Following the arguments from the previous section one can argue that 
the production of identified strange particles in heavy-ion collisions at RHIC 
energies occurs by two different mechanisms in the core and in the mantle 
of the interaction region. In studies with a single thermal source different 
degree of strangeness equilibration as function of centrality is seen ^ |2] . In 
the core-mantle model we choose a different approach. Since the physics of 
the thermal core is dominated by the phase transition we assume a complete 
equilibration of strangeness in the thermal source, this is consistent with the 
data on particle ratios for the most central collisions. 

Unlike the charged particle density which scales with the number of 
wounded nucleons in d + Au collisions, the ratio K/ir depends on the cen- 
trality. Therefore, we cannot assume that in the mantle the particle ratios 
are similar as in p + p collisions. Preliminary data from the PHENIX col- 
laboration ^3] show an increase of the Kj ir ratio with the mean number of 
binary collisions v = A / Vft '"y 2 . This dependence can be fitted for the K~ /ir~ 
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Fig. 2. K~/ir~ ratio as a function of centrality. The data are from the PHENIX 
collaboration (141 IT3| . The solid line denotes the sum of the emission from the 
thermal core and of a contribution from p + p collisions scaled with the number of 
collisions in the mantle (equation EJ . The dashed line represents the results of a 
calculation with a thermal source and a scaled p+p contribution for the production 
in the mantle. 



ratio by the formula 

(y) = 0.062 + 0.025 xi/. (7) 



K~ 

7T~ 



dAu 



Within the core-mantle model the number of produced particles can be 
written as 

AT-m 1\TC 

AT _ ^part e l ^part th 

- ~^—n K H — n K 

N m * N c + 
N* = -f±K l + -f±nf , (8) 

where n el and n th represent the production of particles per participant pair 
in the mantle and in the core respectively. Using n^ 1 ~ (1 + a)n^, i.e. the 



8 



pion multiplicity follows the charged particles scaling, we have 

N K KarS + ^ l + a ) N parS 



Nir N™ rt + (l + a)N- 



(9) 



part 



The K~ /ir~ ratio in the thermal core should reproduce the experimental 

n th _ 

data for the most central collisions. We take — f^— = 0.16 from a statistical 

n _ 

TV 

hadronization model [EH with parameters fitted to particle ratios in central 
collisions. The ratio in the mantle is taken according to the formula 

n el fa 

0, where v is the average number of collisions in the mantle part of the 
interaction region of a Au + Au collision. The results reproduce well the 
experimental data ^1] on K~ /tc~ ratio as a function of centrality (Figure[2J). 
The increase of strange particle production comes from the increase of the 
number of collisions in the mantle and from the increase of the contribution 
of the thermal core to the total emission. Neglecting the dependence of 

the ratio -# or ( ^-^ on the number of collisions we cannot reproduce 

K l" / dAu 

the centrality dependence in d + Au and Au + Au interactions (dashed 
line in Figure |5J). We note that the smooth scaling of the production from 
peripheral d + Au up to central Au + Au collisions (solid line in figure |2J) 
cannot be extended to the data for K + /tt + . 



4. Jet absorption 

Jet quenching has been advocated as a tool for the tomography of the 
dense phase in the collision [16] . Experimental data show a significant sup- 
pression of jets in Au + Au collisions but not in d + Au interactions, suggest- 
ing that the reduction of the jet rate is due to the dense medium created 
in heavy-ion collisions ^1 EH 120! ■ These results were supported by the 
observation of the disappearance of back to back azimuthal jet correlations 
in central Au + Au collisions 19 . The origin of the jet suppression is the 
energy loss of jets or the absorption of jets in the dense matter. In the 
following we use a model of jet absorption j^J which can reproduce the jet 
suppression in the high p± range. 

The survival probability 

/ = e~ kI (10) 

of a jet originating from a binary collision at (x, y) and traveling in the 
direction (n x ,n y ) depends on the density of the matter along its trajectory 



dlp(x + ln x ,y + ln y ) , 



(11) 
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Fig. 3. Raa for charged hadrons and for 7r° integrated for p± > 4.5GeV as a 
function of centrality |22| . The absorption only in the thermalized source is denoted 
by the dotted line. Absorption proportional to the wounded nucleon density is 
represented by the dashed line, the same is for the solid line but with an absorption 
increased by 65% in the thermal core. 



modifications due to the longitudinal expansion and the increase of the ab- 
sorption with / largely cancel out and are not taken into account. Taking 
for the density the wounded nucleon density p = d 2 N part /ds 2 , the jet sup- 
pression at high transverse momentum can be reproduced |21j . Using for 
the formation length of the jet r = 0.2fm and for the absorption coefficient 
k = 0.12fm, the PHENIX data on jet suppression [22] can be reproduced 
(dashed line in figure |3J). 

In a different scenario, it can be assumed that the absorption occurs 
mainly within the dense thermal core of the collision. To illustrate the 
effect we use the extreme limit with a large absorption within the deconfined 
thermal core and no absorption in the mantle. This corresponds to changing 
p to 9{d 2 N part /ds 2 — d cu t) in equation and taking k = 0.8fm _1 in 
equation ()1U|) . We use a larger value of r = 0.8fm, which corresponds in 
this scenario to the formation time of the thermalized opaque core. The 
jet suppression in the most central collisions is reproduced (dotted line in 
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figure E3) • In peripheral collisions the thermal core disappears and the jet 
suppression ratio Raa goes to one, unlike in the data. If we believe that the 
dense phase is not created in peripheral collisions, this means that the jet 
suppression must occur in the lower density mantle part of the interaction 
region as well. This observation confirms the theoretical calculation of the 
jet suppression which increases smoothly with the density, irrespectively of 
the nature of the absorbing medium |23| . 

Finally we perform a calculation taking into account the increased den- 
sity within the thermalized core. Following the increase of the particle 
production we take 

(PN part ( d 2 N part \ ^ d 2 N part ( d 2 N part \ 

9 ~ ~^ 2 ~° [ dcut " + a) ^ 2 ~ 6 [~^~ ~ dcut ) (12) 

Using the absorption coefficient k = 0.09fm (r = 0.2fm) we reproduce the 
data in the whole centrality range. This result shows that the existence of a 
dense core within the interaction region is compatible with the data on the 
jet suppression, if the absorption is proportional to the (increased) density 
in the core. 

5. Conclusions 

We study a model based on the separation of the interaction region in 
Au + Au collisions in two distinct parts : a thermalized core in the most 
dense region and an outer mantle which is not thermalized. The inner core is 
rapidly thermalized when passing the region of the critical temperature and 
above. This includes both the kinetic and chemical equilibration. The decay 
of the thermal core leads to particle production per participant pair larger 
by 65% than in p + p interactions at yfs = 200GeV. The charged particle 
production in the mantle can be taken as in a p + p collision scaled by 
the number of participant pairs. The thermalized core appears at around 
Npart = 40, we expect modifications of the density of produced particle 
and increased fluctuations in this range of centralities due to fluctuations 
in the size of the thermal source; these are not taken into account in this 
study. K~ /ir~ ratio as a function of centrality can be reproduced as a sum 
of a contribution of a chemically equilibrated thermal source and of the 
production from the mantle. In the mantle an increased strangeness per 
participant must be taken into account, following the d + Au data. High p± 
jet suppression is studied in a jet absorption model. In order to reproduce 
the results in peripheral collisions, absorption must occur in the dense core 
as well as in the mantle. The increase of the density in the thermal core 
resulting in an increased absorption of jets is compatible with the data. 
Applying the model to Cu + Cu collisions at yfs = 200GeV we find that the 
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charged particle production depends on the number of participants similarly 
as in Au + Au collisions, up to N part = 120. However, we expect that in 
this case the physics could be modified by fluctuations of the relative size 
of the core and of the mantle. 
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